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dhile the Linear theury hes been the beckbone of ballistic alssile 
design and has served with excelient success in its various foras over the 
P years, recent spectacular aisaiie flight failures have appeared unexpectedly 
end are wmaccounted for by this lineer theory, It is these flight failures 
which now require attention end consideration of the nov-linesrities in the 
fluid foree abd noaent system which costritute to aissile flight performance 
Portia eet eae aE ret cae Bate Sete hate . and dynamic stability. Ges See 
Two types of non-linear flight instabii..y uzve been isolated and 


ddentified eat 
te Son-Lineer Magnus instability, ead 
' 2. Cotastrophic Yewe 
10 fe these tuo Slight instabilities veich will be discussed in the sections 
whieh foliow. Approxiaate aathematica) eodele vill be suggested for the 
evaluation of missile dynamic stability. Both experimental and analytical 


. appsicetions of these models will be sade and discussed. = Ps 
. | . Seek: : . 
“piginelly contained in A paper srepared for Dr. Max M. “unn, Catholic 
" Ualversity of America. 


PM selentific Advisor for Astrcnsutics, Bureau of . rinaice 
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" aaibited flight performance feilures not predicted by the Unified Lineer 


to large lowering tolerances as exeaplified by excessive gun or aortar 

bore elearence. The effect OF these initial launching coniitions vas to ‘ 
Dreduce let, amen values of the complox angle af attack. Jt was, therefore, 

euapected that the aisaile flight $if* culties were dus 0 none Linear ties 

ip the fluid Suree and BONE. ayse” “ON complex angle of attack. The 

eatual alsaile Slight failures theaselves fell +>" wou aistinet grows, 

ooe characterized by aisailes neving large rolling welocity and the other 

by alsstles naving esal) Polling welocity. these tue GFoupes will te 

treated separately in the sections which follow since quite different and 

Telatively independent non-Linearitios are iavolved, a 
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SOU-LINEA MAGI INSTABILITY 


wase certain subsonic fin~stebilised rockets baving high spin vere 
fired cross-wind fros e aoving ship, large and persistent Vobbling motions — 
opeurred. . It vas Pee an examination of the aspaek eriterion for 
dynamic stability, that « non-linear ee shod with complex angie of 


attack aight be the culprit. dsting upon this conjecture a ainor config- 


2-33 
wretions] chenge vas introduced to reduce the Mages somat. Subsequent 


track firing results indicated that the flight instability wea reduced to 


toserable Liaite end, thus, tais critical veapos ves released to the fleet 


few fall service use. 

Thies experience had, ‘ueieen uncovered @ possible sca~-linear Megnus 
peers which Pequires Wi? «fy and understending. accordingly, au 
extensive vind tunee] progran was dadertaker *s ties Nawal Urdnance Laboratory, 
@ non-linear computation program wee undertaken by the Ravel rroving Ground, 
end om integrated flight testing end evaluation progres ves undertaken at 
the Naval urdnance Test station. 

A eiaple empirical aethod for predicting this noe-linear Magnus iaste- 
RAlity wus suggested by the writer end wes eaployed os « guide for eetually 
e@liainating a amber of aissile flignt failures of this type. (ne of the 
purposes of this paper, therefore, is to set forth this aathod and to 
suggest 0 mothenmatical model for substertieting it. 


queai cisear Tieory 


Tligtt instedilities in oallistic missiles are characterized by an 


elaost pure circular pitering and yawing moticn which grows to large values. 
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. In the Lineker case this action is casily traved te an instability is 

eitner the autational or preeesaionel components. ~It 12 the nce-linear case 
whiet new requires attention. As 6 starting poiat it seems reasonable to 
aeeuse “a pricr:* wat tho aissile has om alacei pure circular pitehing 
and yowing action. The dynaaic stability of the missile may, thea, be 
emmmined vith reference to various sises of this special action. (Fig. 1) 
It should be euphesized that this happy state of affairs of alecst constact 

_ magpitede of the complex angle uf atteck exists in rolling fin-stabilised 
Rissiles and shells but does not exist in the case of the pure pitchirg or 
pure yeuing avtion of noe-roliing atreraft. 1 aay be for this reason that 
this approesh has net been considered previously in aissile flight. 

the first use of this assumption vill be 1a specifying the serobsllistic 

force and aomsat system. Al) the Sorces due te angle of atteck will be 

__ spetified by a constant force plus a perturtetion. (rie .} Two forces 


and their moments are sllowsd tc be von-limear. These forces are the normal f 


force dub to angle of attack and the Nagnue farce dus to angle of attack 
‘ amd rolling velocity. Al] other serodynaaic forces and scaents are sssumed 
fo be Lineer in teeir respective variables. decordingly, tne norael fluid 
force and the sormal fluid saoasat acting on 9 missile executing sa alacst 
pure circular piteAing and yeving action say be. written es 
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Roteticcal and wirror symmetry is ssevmed. 


a A ; 

*.* 2y ripe, My? p My ~4 My 

2, ra, ta pe, A? phe 74 M, 

77 tpt, tek, Ms + Me - 4 tM, (3) 

2; 2 2, tapes, Moms Pp Meg - 4 My 

Ro Thea te a Hr ? M TA pelts 
lsWee Wo> ts) 


Tee basic di.terential equations of action written in terms o: 


deroballistic Anes (Fig. 3) are given by*® 


(#,w +d, w) 5 m1 (d-4ug) : ($) 
‘mM, we Mt ww. ary 7° - - —o i 
mf ial: as Be tM uw =a rf 4 foo (5) 


solving aq. (5) Koad 4 ’ iat tor iz acd substituting into 
eae 4°) yields a single az:serentisi equation in coaplex @ which yy 
substitution of tq. (7), ws W rar) e sie and its derivatives 


" pay be reduced to a di..erent a! “equet.on .n real gas 
ce - , 
Seu wW enw = fi, ) 
Sq. (ce) any be uritten with real coefsicicnts ee 


ar ap (%) 
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.provided that the condition serceed by =q- (10) 0 satistied. 
led, - Rs joke 70 129) 
wd N : 


®hor ‘convenience ve shall le ve shall CT 
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where 


my * Revit, o(- Fe - YM) e (-pbpe "26-9 Bop Hes up My) 
(12) 


ie» ches Pree nso Peels) 
+(6 - Fy 6 + Af be “Paes plant - M6~ “Pew “fed ) 


(32) 
m= Rs, hy «(- tiers +P. ‘te tov otto o. e~ rad. "hae ae a 
im Tt 
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The condition «f Ea. (33) may be satisfied either ir aa, =O or ir 
the quantity in the bracuet is enual te sero (1.g. af <i 20 and it 
; i \ * 0). considering the letter situation Sha writing the expreseion 
for Bi, as, 


RN, + - [Sel -"e, tex, * es) 6 ("ere ses doo -"yfe ete] 
In 


(14) 
cae (15) 
ty Brees + eles, ial ad ere PP ty ig) 
we OS ate. 7 ar 
— =+ 7s ar) 
c reFy 
as i oe \ tT - 
3 uM, 


18 my be noted that m, = 0, in@ar,, or if @*7,, 
Welting the expression for um, a, 


—_ 


ange] oC AEE odat ote) (Peo ABE Abe hea] 


= & 
(29) 
and assuming that : 
Bw * aw Ms My (3} 
ylelds ; ‘ 
RN, wz (6-2,)(6-2,) : (22) 
Therefore, it may be nuted irom igs. (15) and (a2) that the condition 
imposed by tn. (3) may be ss°'sfind Af 
eo ? a, : oe oe = e (a) 


provided that the Magnus vorse aod the static monent ere lineer, fq. (22). 
o@ @ay now return to the basic QS: Terentiel equation of e@otion, cq. (9). 
Assuming that e =O, the general solution is given by 


(<3) 


(2a) 


Gelere considering in detail) the solution gives by Eq. (<3), it is 
Gnstrable to obtain expressions for §,, 2 amt ki, 2, 5 1m terme of the 
mMabdility derivatives. decordingly, 4, z 28 given by, 


z, eR) A. lie, a Hew 
Pe “(% Tees lee 6 + i tte aE 
. 4 5 pre \* - Zz ~ & Mf, ‘\6 
t frs-ax6-avs FE) +e - fer - 4 Sev) 6 
A “ YE Mon en 
EE yr Me)- Ry OSE eat, alter, we Me) 
: Tais expression aay be iurther ‘siepliiied if the terms under the radical tz) 


: f6 divised isto two groups, one coateining terns whic are large and the 
ether containing the smell terns. Jalroducing the approximation 


ua Sarg 2th Bek are ea cate eaten Me 
——— ; 3 
‘ : yi Ste ie rs 
Ey. ) becomes, : 


z ad Mea M. ae a Mea 
= ae (4 zt) t + Sar ety = eae a 


ree 
. Boe Ww 2 Mons 
pt ss ye, ) 1. Pee , PSone p Pipe ap “pay Sig 
+ (EE dpe (rer) azr’ an art ary, (9) : 
¢,. = Mik a Wye : 70) 


where the condition imposed by Eq. (1.) Kas been satis.ied. shile Ges 
originally allowed to be cumplex, 1% may be noted froz <q. (29) that it 


/ Se real. To sete of values for), , say be obtained depending upon 
. whieh condition in iq. (ug ) ie satielied, es 
ir ben, um ; 
: » Be (yt) + bs Lala (+t) + afer Qu) 


© fy (ot ry oyete (f-€)- Ar oP ery etew. ‘se sete 
er, ir @e P20 or a 
4, ¥ ae eae “hye ter. alee othe, ca 


sae 433) - 


A © hese piety = (3a; 


—-2 
ra expression tor Ay ia teres of the stability derivatives, when 
bq. G.) is setioried, is gives vy, 


oM, 
Bese) pet Melg ess re | 


W macy W 


rte preted, ite)t “bez Ay 
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: Ora, K. sa, 

3t should be noted that the dencainstor o: the brecke! quantity, (35) 
A. ie identicaily the semiiier lineer theory expression sor the 
autetion ead pre:eseion damping rates where the tangent slopes are used. 

The cumerator om the other Rend while identical in fore uses the secant ° 


slopes for “me non-lineer quantities. (big. 2) 
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a* the eutect, is Speciiying tae : luid Joscen, hag. (3) and (2 ), it 
wee Assumed that the aissile hed alacst pure efreuler pitching and yaving 
Motion. The selection of Luitis: conditions is Geteraining ky and kz 
Steuld, therefore, be such #a to satisry at least daltially this requireasnt. 
Teas, we assume thet, 

ww, Ws W320 
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The eapression sor A), 2 ig. (29, reduces to | 


ia considering the stapicity or » Matile, ve are latereated ja tre 


. CORditicna unter wich the Missile das only marginal dynsaic stability. 


tbat fa, we wah to find @ beand.welow o- AD0e* wh. the ai selle da either 

dynamicaity ateble ur unstedle. Tous we ars ablezested in studying iq. 27>) 
‘oF sonditicnn where \, or Ayis Spprosching sero. It may be noted frog 

a8 (31034) that, usualiy, when one A ts emall tne other Se large. 
clcPdingiy te. | 3t} aay be written se, 


ay - 
tal a ede “s) 2 a, eo 
ms » a 
ae See SS 
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wubstitutang these expPessiona end Eq. (4) tack into fq. @3} yielus, 
me ree 
jar le “she + (We «,) 


¢ ig Pe 


te 


(36) 


where asd, tam 4, wo 
or duh, am 4 0 
. Since F, 18 asonctated vita A, , and rz with A, , ve aay aiaply veite 
our solution for the antics as, 


g v: fee Sea, aed 
la} «(ur+ewye +W, { ae 
é 


2 , eo 
(a) 
where " 4 
We Cl we (42) 


Tals ie, therefore, the besic equation for the alacet circuler motion 


of @ alssile ected upon ty B non-linear Magnus eoment. In the following 


Section this equation will be enanined for evaluating the dynamic 
of the atesile. 
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Distussion of the Quasi-Linear Theory for Magus Instability ~ 


Our perpsve here is t+ dissuse Eq. (41) with a view towards deter- 
mining both the conditions for aisslle dyuanic stability aad a measure of 
a cee oe dynamic stability vaieh « aissile hes for given flight 
eenditions. 

‘The conditica for dynamic stability is this aun ais gee by, 

o> ne - [270 +%) pret se) + a (43) 
Thus, a0 6 starting point for stutying iq. (41), 4% seems reasceable to 
begin by exanining the case of negative A. First, it may be noted thet 
warn \W=0, then /ti/-Oend alos wren y=) , /u7/-* 0 » Both of these 
Pouused cases are eoventially **= familiar linear one when only one arm existe. 

Of baste interest in tats mcn-Linear solution, tq. /..), is the fact 
that © trin-like condition oxiets mush 19 represented by Wr. ta 
wonsidering alssile stability, it is i-portart to determine the size of this 
peoudo tris as compared to the basic motion, W , about voich the pertur- 
bation, Uf, takes place. Thus, values for “iy are given by Eq. (42) 
and plotted in Fig. é. ‘ 

In the familiar Linear case ween \<O dynamic stability was indicated. 
Bovever, in this noo-lineer sclution it may be moted that if Wy > W esd 

d <0 , then ao divergesce from the basic motion results thereby 

Indioating dynante tastabitity. Cf course, if W, <W ant ) <0 


dyeandc stability results 


Tae wee Of the terme Sdynante otability®, ete. in thie now-Linses, i. 
Case should be defined. 4 hesic action Pepresented by W is perturbed 
Ny the eidition of Of | Fig. (8). IF the remalting action, [ti], tends 
to reduce in sise approaching Wor « value less then W » then the 
motion 16 termed “dynamically stable” to: that value of W. 

Tee various possibilities for etability end jnatebi lity in this none 
Uneer case, iq. (42 ), may be perceived by considering the two besic 
caves vhich are Fepremeated by Wr cw, Fig. @, ant WW, He. é- 
The corner 1@ obtained when Mh 48 positive ant the letter unen’ Mh ia 
Begetive, (see ta. ( 42) and rig.g: ) 

Two possibliities for positive . 


A exlat. eben A de negative and 
{2 negative, dynaaic instatility results. le Fig. § these twa 


possibilities are pee Nes os Me “dpae posetbliity (2 <0 ang ie “), f < 


the transient part, (  ) - 
towerds \t 


» SRFinke, and trv 2. cote) aotion tends 


+S WW, thus ratieetiag Gynamic stability. In the second 


possibility (A>0 gag 27> ) the trensient part ova, and thus the 


Rotel action tends to increase avay cron 4, wd: W , indicating dynamic 
instability. Ppa 
. e 
= possibilities :or negative * , alac exist. shen * i hegati ve 


: ae positive, the afesile is dynamically unstable, Hovever, when 


‘ 1 positive end * SB negative, dynamic stability resulte. This 


satter possibility is o: particular interest siece it certainly would not 


be Goticspeted from any linear emalyaia. 


and 


These two possibilities ‘are 


Z ¢ 
Pepresented in Hig. 6. In the fdvat NO em NO), toe transient 


bart shrirke in sise end this the total, action grows avay sree ' and 


ao 


cece ee mage nage scoters tte olay Semana sn atetaninatan «tae st = at ete tne mae en rt 


towards W.>W , indicating instability. Is the second ()>O and 
d* CO ), tae transient part grove in nize aad thus the total action 
serike toward «OW , indicating stability. 

aan \'s 0 and ‘ie mgative, the pitching and paving motion will 
be stable in the criginal pure circular aede neither damping nor expending 
since W, =W. The stable asdes of etrevler ae oe *Lamdt Cycles*) 
have been previously inve-tigeted by other writers. 

mes d.*0 and A ie positive the aetion will grew and thus be 
wmatebls. 
: Walle negative values for */y nave act been specifically discussed 

oe ililvatrated, all geceral stateceata still apply. Por ecapleteness, 

ome typical esse ie represented in Fig. 7 & 8 and is now discussed. Wheo 
XO ad vo ond NTS! 3! team Wy OD aut te trenstent motion reduces 
In sine; the total astion tends toverd Wo asi o/s W ant t9, tous, 
i dynamically stadle. 
Ie quamary it hes been shown thet the missile is dynentcally stanle 


vees perturbed froe the basic action for tuo casse ven 
<0 und 


A 

mm 60 od co. 

» and the aissile is dynamically unstable in three cases when 
d 


aw 

e 
a] «oO amt A 4D 
. eo 
a » %O ad AX »o 


aM 


> 


0 and A =0 
e 


& Liait Cycle exists wnen \ +5 


md gco. 


Clearly tas sign of ) 10 ect « asesure of missile dynamic stability 
im thie nemLincer cose as it was in the liscer case. In reviewing 
Bes. (ab) = ( K9) 18 10 feud that the aim of y" io, hevever, a 
eviterion fer dynenie stability. Tams ve Sind, vith seme easement, 
that the reqyuirenant fer Dynemie Stability is this sen-lincer ease is 
es 2 


<o 


ae Lae | (Mp+u M,) uMpy 
4 ae Fe de ce 


‘Tae condition ouggecte tant a plat of yn wet eW an 
imediately reveal an evaluation of the Qymenie stability ef the aisaile. 
Tale to thoenpirical axthed vaich ves originally suggested end vaich the 
Presceding anthenaticel andel comns to cubstentiste. fy evelusting )” 

fer cach value of W using tho setae) sm-Lineer soretynanis coefficieate, 
tae plot is made. 18 sheald be recalled that tue exproseiens for ’ 

euiet depending on vaetner 2, of *7,, Thue, two plete sbeuld be 


ante. In heeping vith elessical usage the fester eireular motion vil) be 
called “nutation® and the slover, “proceesion’. sncoréingiy, regieas 
ef ahesile non-linear ¢ynante etebility any be revealed ia such plete 

GO given Wy Pig. 9. 


Set enly any the dyncaia stehility be determined bet, alse, the 
Muamah of Gyrende elahility any be evalusted by using the non-linear 
theery. In the stable case, by eubetiveting |i] = Male wow 
inte Bq. (43) yields am enprescios for Notion Half-Life as 2 fnstion 


os Re; aon z oli tL Z | G2) 
ue” dy, 2(w+ ew) ; 
Aue ¢ . w) _ 
v ‘ ) ? 4 
x= bii- 
ddue | Phu n cow ed 


iis sie natu sans: os agemeaton’ Pe WAdin BeAdOALIEN eng ba 
evtained ty eutstitutag fix] * 2 [6%] tote fq. (41) 09 


a es ee 
2. =v 4 ae tt Ww 4) 
Pd hue ee ap 
w= L| l- | 
d dap Bane dence ($5) 


Aa @ reealt of the preceding enalysio it bas bean fvund that: 
(2) tee condition fer dynenis stability is X(w)<o 
(2) the ensunt of dynante stability ts gives ty the enpreseime 
fer Motion Half-Life end Deuble-Life. 


dgyplications of the Quaai Liscer Theer7 


applications of the Quasi Lincar Theery for Kagues Instability have 
been aade te certain atesiles vaish, in Plight, exhibited « pecsible Neguus 
AmatebiLity. The appreash cuployed une, firet, te make plote of Nae lW) 
fer the original design, for asdified designe, and for nedified lemebing 
conditions ant, thee, enlect these conditions for vaieh pocitive valuss 
fer ie wore entirely aveiéed. Subssquent designe based en tats 
procedure vere all ebeorved te be éyneateally stable ia flight. ade 
the epectesuler cuscees of these Limited enperimmatal applications ef the 
Qheery indicate prenise, the full extent of ite weefulaces and of ite \ 
Liadtetions te cocentiolly winrws. Sines ve are dealing uite o highly 7 
ten-lincer problem, great ooyic!~n end care cbeuld alvaye be enercized. Ss 
In atdition to the experiasatal epplisetions. orresieen of the , ; 
OUDILity prodiations of tals theory vere ante vith resalte free smericel 
tategretions of the cnset nen-Liseer differentia! equations of ution on 
{ho Bevel Ordnance Ascserch Cosputer ot the finval Proving Ground, Dehlgres, 
Virgiale. Thess cenyerisens of the predictions of the Queal Linser Theosy 
with the BOR resnite are given ty Dr. C. 3. Coben and 8. C. Inbheré in 
aga wise Tuas> task ‘Soiadas aus oacaie Ga Saeed 
Piaally 10 ebsalé be capheciond that the Quesi-Linser Theery ie vet 
" fatended to predict the gueeral detailed free flight aptien of atssiles, 
particularly vhes tet® sutation and precession erme are present. Rather, 
the tasery euggecte o anth~4 for quiekly and cocily evaluating nissile dynente 
wrediiity when ested upon by 0 sup-liscer Mageus anacst. Por initial design 


evalantions and 0s 9 guide fer “quick fizee*, it 19 strongly recomended; 
bewever, 18 19 not intended to replace fical sumvical integrations of 
the emast differential equations of aotiin os « high speed ccaputer. 
Muperionse, tame fur, has indicated that the best results end e full 
pavsioal wederstending are obtained by using both. 


Introduction 


A unique flight instability ast accounted fer by the Unified Liseer 
Tasery erty Soe Greet Linmer: Teeeey’ eae’ eesestenally ‘shane wet tn the 
: flights of « cortala fle-stebilised alesile. This particular dyeanic 
. '" tpgtabdLity, ween it eoourred, vee first eberesterized by the failure ef 
the aisaile to plok wp its full steady state Yelling veloetty and, then, 
by @ catastrophic growth of the pitching and yeving action. Initially 
the relling weloeity tacrenses due te the fin cent. Hevever, vhen it 
Feaghee @ value equal te the sutetios frequeasy it helds eenstant and “loeke- 
ia® at tact particular value iurtead of seeking the mesh larger steady state 
it value. (Fig.10) Initially the olse of the pitchise en? ,awing aation reduces. 
Novever, when the rolling mations locke-in, it beging to grew end may very 
eoca reach extreme values such thal the missile is flight lboke mere like 
; @ propeller thas an arrow. (Pig. 11) This waique flat spin hes bees terned 

“coractrepaic tau®. 

In attempting to find @ physical explanation and possibly @ theory for 
the prediction and future avoidance of the Catastrophic Yew, the rolling 
motion will be considered first. 


Rolling Motica 


Tae rolling action of missiles hes been classically considered the 
sinplest aode. Its hnovledge end control, however, to essential in the 
@aided aiseile in order to achieve guidance and in the ballistic missile 


19 


ia evder to ashieve flight stability and acewary. Por aiseiles firing 
ot anal! anglee of attack, the Linear theory hes perforned nai 
la evder te critically investigate the possible existenee and natere of 
@hnereal rolling actions, pratisulariy of tbe type eseorved in catastrophic 
flight, « simple Sie-stabilized aede!l (Basic Finser) wae sting-mevated 
Svan: sired. th he ‘issiccah Sirsne' 6 Disnhendé-ciad Simesks Ths 
veriatios of the steady state relling welovity vith mgle of attack for 
angles fren U* te 3” vas msewred. Typical date are given ts Pig. 12 
were it ie meted that there are large and alasst distontinuous ehenges 
ia the steady state relling velocity with angle of attack. Ales, it is 
noted thet the missile any rell in either direction and, if the aisetle 
ip stepped, it may remain stopped. , 
|" oe ecnventence tne tu. mi abnormal relling astions cheerved during 
Sisise aapherthaty’<oihG aie Uukahion Seah tassel eck 
Slace, over o certain range ef angles ef attesk, the aiseils can perfors 
-elther setion, it wes suspested that the tus phonemes are basizally 
independent and muy be treated seperately. fer our present purposes, it is 
tae failure of the aissile to rell (roll lesk-ia) vaieh required ealysis 
end eaplanatica. a 2 

Tue bey fosters characterising roll lesh-in are gleaned fron the 
euperizeateal data. Firet, there ie a definite depsndence cn beth the sngle 
of attack and tae fis cont angle. Belew e critical valze of angle of 
attack, the aede) will rol.; abewe that waluv, the aedel may fail te reli. 
In the range ef critics! values for the eagle of atteck, there is a 
eritieal value cf the fin cant angle below uaieh the andel vill fail te 
roll asd above waich the model wil. roll. Sesendly, when the atseile fails 


2 


te rell, the rell orientation of the aedel vith reapest te the plane of 
the angle of attest, ¥ , verics vith the angle of attesk. At the critical 
angle of attesk, the rel) ericata:ies of the model is in the neighborhood 
of 67 1/2. ta the angle of ettask 1s inereased the rol) erieatiticn 
tends toverd 45° on « linit. 

These tue eharecteristics of rell losk-in suggest the possible 
existence of a rell senent waleh depends on angle cf attack and aissile 
Fell orieatetion. A review of the litersture reveals that ousk 4 see-Lincar 
sell aansat is know to exist oan< bes bow srccured quite early oa airships, 
torpadoos, and alreraft, “for o cresifere fined aiseiln, this *Indueed 
Rol) Monsat® tahoe the general form 


Lita) = Cot Cam tyeav'Sd 


(56) 
Ga Ue 
ant to plotted ia Pig. 1. 
"= la are interested is deternining if the addition of this sen-Lisser 
yell saneat te the Linear rell equation cam aocowst for the phenomena of 
roll leek-ia. Ascordingly, the equation for the rolling setion aay be 
rewritten as ; ‘ 
L(Se)+L(p) +L (ya) = I, p (57) 
ware 
LCS) a Manan: due te cant = Cs, (a) 5, 40d (58) 


Le CP) ® Rall Damping Hemet 9 Cy, (a.m) ([BL]$pv'Sd (59) 


L(Yd): = Moment 9 Cy fet, t) a Lam] apy ad 


In the iad tee) cose of roll lnck-ia pro, GeO , and 
of = const. Thue Sq. (57) reduces te 
L (Sa) + L(t) 20 (a) 
In any particular wind tumel test of roll loek-in, 


d, = tongr ol = tomir 

Tens te entity Bq. (61 ) 1% is only necessary to fied ¢ velus 

of Y fer vaies P 
<5 an" [- a | (62) 
Via! 

The tures possible nee-trivisl casse are illustrated in Pig. (14. 

In Cone 1 the roll nement due te cant ie greater than the induced 
reli scmeat for ell veluse of the re:t orientation. Thus, there ie no 
welus of oc satisfying Bq. 142 ) and therefore te. exte! cannct lock-in 
eed vill roll. 

"Im Cose-Zthe tue roll moments are equal and opposite at Y-# 67 1/2°, 
and thus Bq. (62) 1d satisfied and the model ces lock-in. Thie 1s seen 
te be the “critical angle of attack” for thet perticular roll momeat due. 
te cant. If the angle of atteck is slightly decreased it le acted that 
the sedel cannot lock-in aiace ve have Case 1 egeia. If the angle of 
attech ie ineressed then ve have ancther example of Case }. 

le Case 3 there are tue values of X for waich fq. (62) may be 
satisfied. Por the value 67 1/2" < ¥ © . 90° otatés roll instability 
eniste since o slight increses in \{ yields « roll aemect due to cont 
greater thes the induced rell someat which tends to couse a greeter 
increase in and thus further static mbslance. If ¥ te slightly 


s 


** 


Pedwed © siniler wmotable situation existe. Mevever, valuce of 

we ¥ $ oe ees te essay eee me 
Pel) lesk~in ie eapocted ia this range, ° 
, si giv Gasiau ca a lentes ee 
voll tein engle are thus given by fq. (43) by simply inserting the 
payelea) peren:*ere invelved, 


Lee, ac “a Be (63) 
tye 
Ca ba F 
Vien * = am - ee 45°S Nig £6Th 
(44) 


vises Ca. ond Cy | aay both be nemelinear tn << nd Y. 


10 wae aloo observed in the vin’ tamesl *>..0 tnet the model vace 
canited ted nal) coe! listions in rell abewt the reli trian angle. the 
froqeeney of thee cosiliations any be predicted fren iq. (67 ) ty asountag 
@ Mesar variation ta (1a) for the anall ¥ renge involved. 


" geeordingly fq. (37) any be written a6 


bir ey a s 
Stu ve Van, ” 
a N, -¢ [4,18] tev'sd ] —* 
N, aa ae & tpv'sd] aa pe (67) 
Ls + Blogs tav'sa] cat 


a 


- bel TS=3 = Pa 
A, Se 

d, ate Lm -4H, * Att @, 

ge s = = 

ws *, ~ N 

K-EAUwEY eee 

a. * —y" 5M) &j Zz 


Tee frequency of the roll cesillations is gives by 
O48 (Ea tevsa)™ 
and the damping by 

$= se idplb}tev'ss) 

‘‘siathhinn, particularly ot very large angles of attack, the aedel 
was cbeerved to cecillate in's divergent fasm:.c a’ beets to roll ty 
iteelf. Thie type of reli performance is a'=: a...cmted for by Eq. ( @) 
if the roll damping sement ie positive or wademping. fhe above analysis 


of the experineatelly cbeerved divergeat rel) ceciliations suggests that 


the rell damping soment aay po non-lincer ia the angle of attack and the 
Felling velosity of the atcsile. 4 sce-Linearity of thie type aight vel 
account for the roll speed-up phencaena meutioned at the outset but act 
specifically otatiod tere.” 

Time, the introduction of the noe-Linear induced roll mmeut inte 
the classical Lineer roll theory wnich containe the roll ‘moment due to 
cant end tae roll senmet dus to rolling velocity appears to yield an 


explanation of ead a prediction method for the rell lovk-in chservad 


ia the vind tame] isots. 


. 


Rew that the wing tune) rel] loek-in phenomene appears to be 
accounted for, the free flight missile rel) case aight be profitedly 


considered. This aissil> case presente tue additional complications. 
Piret, the missile ie rolling and ssscad, the siseile's engle of attach 
aay be changing. It is recalled thet in the miscile case the rolling 
velocity failed to increase to ite design value ead it stopped increasing 
at a value ef relling velocity equal te the frequency of the aissile's 
pitching end yering action. Sow when a aissile is relliag at the seme 
rate that 1t to wobbling and tf the wodbilag is in the same sense as the 
Polling then we have the sperial ease where the alssile is net changing 
ite rel erlemtation with reepect te the plane of the complex angle of 
attack or, eald ethervice, Y te a constant. This type of pitching ant 
yoving asticn io called "Lams “tien” peceuee of ite analogy to the 
mation of the meen. ; ; 

If pare luper pitening and yuwing aotion existe, where the angle of 
attack ie o constant, thes the analysio used in the wind tunel case may 
readily be extended to the missile case by simply adding the ell damping 
manent. 

Oy. (’97) ang be written as 
L(8) + Lire) +L(ya) =O 
vere &,, 40, % , ant Y are conctente, the eritical angle of 


atterk io gives by 


fo 


Car 


Ch ba -%, (84) | 
‘Ne am FY 5 


as 


(75) 


and the roll tris angle is gime Wy 


ol , 
Yon #4 [25 (6-508) oom 


In Pig. 14 the elution to iq. { Y¥) is represented. 

The phenonenen of rell leek-is ae observed is the sperial vind 
tamnel toate and in full scale free flight of « particular niseile «appears 
to be traceable te the actions of the men-lisear induced rell manent. vhen i 
thie aameat is introduced isto the elaseical theory fer the pure rellisg a 
getion, a method for predieting ths critical angle ef attest and the roll 
tria angle is evolved. 


Pitching and Yawice Mo**.. 


3t now remains to anaver the second question, “way the catestrophic 
grovth of the pit:ning end yowing metical” It vas eateblisned in the 
previeus sectics end borae out by the flight performance data that viea 
rel] loek-is cecurred, the pitching and yaving sotion is of the “Lumar 
Type". Specifically thea, we are conearned about the catastrophic grovts 
of luer metion. For « statically eteble missile, the Unified Linear 
Theery indicates that lunar aotice any occur im two distinct cases. One, 
when pure sutaticn (1.0. A; ) exists and the autation rate is equal to the 


roll rate; the o.ter, vben pure trim (i.e. “; ) exfete. Is the letter 
‘ease lumer sotion existe for all values of the aissile rolling velocity. - 


Of evaree, a conbination of beth of the couse yielding lunar netics can 
exist when the relling velosity is equul te the sutation rate. Altiough 


tbe conbination aay be mere realistic, we shall first conaider the tue 


peosibilities separately. 

Tue dymanie stability of pure mutational aetion in the liner and 
ia the ace-linear case was treated in the previous seetions. ipplying 
these nothede of predictica to the aiesile in question ecspletely fails 
te evcomnt for the chesrved cetastropaic grovth. The sise of the tris is 
well knows to be 4 fusction of vse rolling velosity of the sisaile. 
Meaxioum asplifications of the nor-reliing trim cccur when the rolling 
welosity is equal to the autation rate. Thie phencasaca is know as 
“hicsniahee testadilig? and initially one would certainly guapect it as 
the cause ef the catastrociuc gr. th °° the pitching and yewing astica. 


_ However, thewe amplification factors can reat'’; ue calculated from the 


linear theory and even apprezisate valuse can be obteined when noe~lisear! tics 
in tae earodynente fovee end aoment systes iekeke These calculations reveel 
thet am explasetion for catastrophic Yew cannot be found in current linear 
OF non-linear theory waich contains the classical static aad dynamic forces 
and moments froe esrodynaaics, hydrodynaaica, end ballistics. <ieerly « 
nev fluid force or forces are required for an understanding end @ possible 
solutian. : . 

@hile dependence os the roll oriestatioa io immediately ruled out ; 
ia tae lineer ceee, ite essential role in providing an eeplanation for roll 
lock-in suggests that it mignt profitably be considered ia the search for 


& force end aopent which might account for cateetrophic yaw. Wiad tunnel 


teste de roves] two ebditians] effects of reli crisatation. First, tas 


"~~. aeeal forse and its manent ore andified as indicated in Pig. 15nd, 


seecad, & side force and mare are found to exist and are indicated in 
Mig. . : 
The nepal feree and ite senent any be writtes* 


Bit) = op w tey's + Cay 7 (am 28) 4 PV'S! (78) 


4, (ar) = Smo w tv'Sd + Ce A we (wm 24) bo wa (79) 


where the first tere in the right aide represents the classicel fors 
waich in iteelf aay be non-linear in angle of attack and where the secosd 
{ere represents the variation dus te rel) orientation. Hers, also, the 
ecefficiont may, (2 ths geaeral nace, oe 6 fuetion of angle of otters 
end roll erleatatica. 

Tae aide ferce end aeasat any be oritten 


Ziyw)tadyy (am F0jGF » Cn wee ti)/$ pV'S (80) 


M (iw) 2 Mea { ame 48) ur: Cmca ae thy boy Sd <81} 


The equa@icns (1.0. Eqs. ( 78 0 79)) are applicable to @ crmeiforn alesile . 
_caly. Sintler expressions any be obtained for other rotations) asymmetries. 


In evder to doternize the centributions of theee nen-lincar ferece enc 
aenante to the flight perfurmaner of aissilee and specifically te the 
Gyrania stability of the special pitening and yoving avtion, lunar action, 
At 9 neccesary to odd then to the classical coredymanies eysten and te 
investigate the nedifie! equations of astion. Clearly a gueral selxtice 
Cf the complete equations ef motion containing these new nea-linear terns ; 
fe net possible. Reverer, if ve confine cur attention to the special ease Ss 
of lumar sstion om appreniaate solution is possible. 
Taking the perturbation approach caployed in the previous section, 
the total force and mumaat, 2q5. (1. 8), aay be extenddd oo 
= a 4 a4 we » < ~ 
Zz =[(3, +8) (3, hw] x rf +l.w sof (#2) 
Mm [OR lw ec +H vw Le ihe ‘Mw eM. 3 (a3) 
oes wl es a4 
where ” ’ 
2 ne ® a (am2t) es Bw onl) {8} j 


a 


A ela (as) 
Aw Fg (ama) te Bw (name | | 
A ‘ 
Cg : Mw (omy) ~4 Mw (se 24} (96) 


yee ub (amt Pap ( nam 21) (e7) 


in the general eees OF and ¥ very vite tine and the eusfficicnts 
mag very with ig and { . Sevever, in order te obtain an appraminete 
“, salution of engincoring value it is nov assumed thet the coefficients are 
constant and thet Y is appreninately constant. secordiaely, Bw ie 
Aye , mt My, ave therefore constant end substitation of Eqs. (82°83 
inte Sqn. (5°G) yleldstne differential equitions.cf action. 
Tee bende form of igo. (5 6) after oubetituting qe. (62. 89 Loa 
Pennine vashanged enc thus tae anthod ef approach, the solutice, the 
analysis, and ten diseussion of the previews section all appiy equally wll : 
_ here.® The now dynamic stability factor, therefore, teccmss 


Set ae 


“ 2 aff f , 
ee Pee aye SEO SS [hae Malet 


(8) 


SP» oc ee Por 


The requiresent for dyneaic stability i. ela:ned fros Iq. 60 ) 
and is 


w <O (f,) 


The contrituticae of the roll depeadeat forces end soments to the 
dynamic stability of Lunar asticn may sew be considered using Iq. (#8). 
Por aiselles with canted fine roll lock-in was saown to occur im the region 


since ¥ was cammned constant it follows that p* @, throwghent this 
developanat. 


SS YX < ora’. rans tas first tore tn hq. (88) ie redored in . 
tise by the rell dependent foree, The effest 19 noted to be destebilising Pe 
in o finestabilioed aissile. The last tere in Eq. (88) ta the anjor 

contributor to stability, and it is meted that the addition of side mesct 


May indeed have a catastrophic effect on the dymeaic atebility of Lunar 
watica. 2 


Disewssion of Catastrophis Yew Evaluatica 


The -aathematical analyses of roll losk-ia and catastrophic yay suggest 
@ aethod for possibly evaluating this type of flight instebility and : 

' avoiding it is missile tesign. The aetmud considers thé alssile to be in 
pure cireuler pitebing ant youing sotion of amplitade W. the first etep 
ie te deternive if rel) loek-ia so wisthie and, if se, te determize the rol) 
Losk~in angle. By fluid caleulat! cas, wen vesible, or by experinental 

Su ade) Vesting in vied tunnels, vater tummele, sercballiatic ranges, ete 
. values tor Cy 4 Cy, ant Cp may be obtained. These values 
ween oubetituted into Eq. (77) yield ‘values for the rell lesk=is angle : 
Ge 8 fastion of the complex angle of attark. 

The next stop is to compute Ne fone Bq. ( 86) ty wotng 
theeretionl or asperimeatal values for the required forese end mensate sot 
and by weing the previcualy doternized valuce for ¥ (W) 

Time a plot of Ne es a function of the size of the aacunsd ctreular 
aetion, W ¢ any be ande in which the regions of possible dynanic 
inetedility will be revealed. Yalme of XO may be avoided (1) ty 
redveing the initial leuneting emt flight conditions oush that waeteble 


- 


2 


waluse of Wf vill set be encountered or (2) by configurational or iooctie 
redesign of the missile to assure, <0 for all anticipated values of W. 
It 4s clear that either the eliairation of roll lock-in by reducing Che 
or the elininetion of the detrimestal side sonent contritutions to \’, 
ack puna Satestropiic Yew. it is noted that roll lesk-is is « necesserr 
wet net a sufficient condition for this instability. a 

Walle 10 fe recognised that « aissile in flight vill sot gmerally Bree 
ewe a pure circular aotica, as asouned bere, one would certainly have 
eensidereble leek of confidence in any alnetio design which was indicated 
dynamically wstexle in this siaple mode, The requirenmt thet the aisaile 
be dynamically stable in tale simple pure cirouler antion eecas nsndatory 
if decigne cf aargina) stability are to be elininated. Thie my be 
particularly true vher 1% is recalled thet all flight dynenic instabilities 
_ obeorved by the writer aave been observed to be uitiaat *y of the pure 


eirculer type. 
The emtenaion of tais approach to sore complicated astions is left 
fer future otady. Applications, thus Sar, of this Cotectrepate Inv Theory 
are presising and indicete tnat this simple, ensy, and quick method should ‘ 
be of anciztanse to the design engineer end the flight dynanicist. 


Applications of Catastrophic Yev Theory 


Linited theoretical and cpietaacts applications of the Catastrophic 
Yas Theory teve bees aade. In the experimmtal case, wm air dropped aissile 
; une considered vateh was cosasioually observed to have extransly large 
pitening and yoving astione erccapeniod by 0 failure te rel) at the designed 
rete, The motion vas diagnosed to he Catastrophic Yaw and, guided by the 


3a 


theery given here wind tunnel ests vere carrisdet.Values for the required 
aerodynamic coefficients vere obtaired for the basic configuration as 

well as three acdifications, Dynamic stability plots (i.e. ‘e (w)) 

were made for each of the four configurations in order to evaluate their 
velative dynamic stability. Tes full scale designe of each type vere 
constructed and flight tested. The detailed dynasic performance of each 
type was chaerved.' The reaults from the observations of flight stability - . 
om eaeh type were in the seme order of evaluation given’ by the predictions 
of the theory. The flight teat prograt is ieported in reference (41). 


MORC Study 
Tn the amalytical case, the exast non-linear differential equetions 
of action containing i-“uced roli momen. 44 *'de moment vere coded for 


nuserical integration on MAC. The prsow 1 the NCRC calculation ves 
twofold; first, 4t was deslred tc compute the flight performance of « 


” particular alssile when subjected to various initial leunching conditions, 


and second, it vee desired to compare the predictions of the epprozinate 
Catastrophic Yaw Theory vite the results from the integration of the 


“enact equations. 


Approzigately 200 complete six degrees of freedca rigid body 
comg :tations vere carried out on SAC. All the in-puts (ecefficients, non- 


-Linearities, fin incidence eles, physical parameters, initial copsiticss, 
5. 


ote.) were systemically varied. Three cases most clearsy reflect the 


eesential contributions of the induced roll mement and the side soneat 


\f 


te alssile flight perforanace ani stability. These three cases a 
illustrated in Fig. 17. In order to provide = saximm opportunity ‘or 
Catastrophic Taw to occur during the coaputations, initial conditions 

_ were selected which represented Lumar action. in order to obtain a 
eonplete biatory of the missile's perforaance for any given cuss, 
ecaputetions vere run both forward and beckvard in tine. ds @ result, 
for these three cases the angle of attack and th. rolling velocity are all 

- identical at «@ specific value of tiae not equal to 0. The “linear case” 
represents the angle of attack Aistory and the polling velocity histcry 
ween only the elessical linear serodynanic system is enpinyed. Here it 
ie noted that the roll performance is regular exhibitirg no tendency to 


lock-in at the Mutation rate, end also exhibiting a pitching and yeving 


, motion that is well daaped “sre senting the dynamically stable design. 
shen the induced toll scmeat is added to the Linear * culation, the 
mations labeled “induced roll moaent” are obtained. The ritening and - 
yewing motion ts observed to be essentially uncharged as would be predicted 
froa the Catastrophic Yau Theory. The pollirg motion, hovever, exhibits 
Poll lock-ia end later roll break-out. Seth of these characteristics 
are predicted ty the Catastrophic Yew Theory. The roll lock-in oecure 

| wren ¥ lee constant. fell break-out oceure at thet value of the angle 
of attach for which the induced rol] amet is no longer larger than, 
or equal to, the rol] moment due tc cant of the fine. Onder these 
circumstances the rol: soment dus to cant Ls able to be effective in 
eumeing the alseile to increase its rolling veloeity and tend toward tre 
required steady state value. This esse appears to bear out the predictions 


a" 


@f the ecatribution ef the addition of the imtused roll acmat es 
afferting the rel] perfornenee of aiseiles. The last case represents 
the angle of attack aad rolling histories uhen both the induced roll _ 
gquasat end the side moment are added te the linear formulation. Hers ve 
observe Catastrophic Yew and complete roll lock-in. (Pig. 1840 ineluded 
te illustrate the tateresting manner in vhich the roll orientation — 
angie, % , changes to tine price to lock-ia.) 

The three cases clearly indicate the dependence of these dulatttine 
in contributing to the overall performance of the missile. Without the 
indwwed roll apaent, rol) lock-in could sot ocewr. Without roll lock-in 
the side masat could mot ast ecastently in the Magnus sense end cause 
wdemping of tne pitching and yzrving motion. In order to evold flight 
instabilities of tais type, °: ts clear that two obvious avenues of 
approach ere available. One approech is to reuse **. cide manent to 


' @ eiae sweb thet it vill sot produce the dymemic instability. The other 


approach is to reduce the induced roll monest ao that roll lock-in does 
mot cecur. 
A third but net so cbvicus epproach vould be to except the eetion 
of the indveed roll moment producing fol) lock-im, ard also aecept the 
side momeat, but obtain a value for roll trim ia the regice 22 1/2°< X,< 45° 


‘ge that the side moment vould thes ert as « strong etebilising influence. 


& value of 1, in tais region may de obteined by elaply reaoving the angle : 
ef cant from the fins. It is essen free the Catastrophic Yew Theory that 


Fell loet-in vould then occur in this regics. Calculations carried out 
om BURC vith end without the roll cast angle are given is "ig. 19 vhere 


3 ’ 


it is asted that although reli lcek-ia cecurs whee the cant angle is 
: 

venoved the dynamic instability due te side moment does not occur in 

the pitehing and yeving wtica. 


‘Comeludiag Remarks 


the results from both the experizeatal prograa and the computational 
wtudy using the exact equations lend considerabis conficence in the 
qualitative predictions of alssile. dynamic stability using the Catastrophic 
Tew Taeory. 

In addition to providing ‘the aerchallistic design engineer vith « 
quick and eney method for predicting tne bounds of alasile dynasic 
etebility, the analysis has stimulated the development of vind tunel ant 
eercballistic range tectuiques ?rr the experimental determination of the 
nevensary eerodynanic coefficients, aad also bas peev'is. general non-lineer 
equations of miseile actice in ai. six degrees of freedoa which now are 
coded for ready computation on MRC and are available for amy proposed 
aissile ei 

The ever incressing coet end coaplezity cf alssile veepons syetese 
Fequires an extensive use of this new capability which ‘dies the precise 
computation of aissile flight performance and thus system evalue’ion long 
before any detatled mechanical or electric design or construction is 


vadertaken. 
dith this approach unpromising weapons syateas designe may be 
eliminated early in the research and development progres and concentrated 


Sg eh en 


taphasis any be plased "a priori’ on the meet promising eystens cad 


‘their wee peodlens. 
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